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Y,  SUWCART 

Until  recently,  the  one-parameter  singular  expression  for  stresses 
near  a  crack-tip  was  widely  thought  to  be  sufficiently  accurate  over  a 
reasonable  region  for  any  geometry  and  loading  condi tiooe.  This  view  has 
been  fast  changing  due  to  the  growing  evidence  that  the  indue*'*)  of  higher 
order  terms  can  significantly  affect  the  solution,  particularly  v  So,  certain 
biaxial  loading  conditions.  In  this  context,  the  present  paper  /xa  mines  the 
strain  energy  density  criterion  for  fracture,  and  the  consequence!  of  the 
assumption  of  a  1/r  energy  singularity  in  the  formulation  on  Its  application.  It 
is  found  that  this  assumption  imposes  a  rather  severe  restriction  on  the  region 
for  which  the  criterion  is  applicable,  and  that  its  application  on  an  arbitrarily 
selected  ‘small’  distance  from  the  crack-tip  (a  procedure  which  has  been 
adopted  by  many  experimentalists  and  finite  element  analysts),  can  lead  to 
erroneous  results.  _ 
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NOTATION 

A,B,C  Fint  three  coefldeats  of  tkc  (trie*  ex^niot  of  tke  strain  tiet|y 

density  faactioa 

a  Halt  crack  length  of  specimen 

h  Halt  width  of  ipccimca 

E  Elastic  modahu 

F  Faactioa  of  elastic  modalos  sad  Poisson's  ratio 

A  Half  height  of  apedmea 

i  y/=l 

J  ./-integral 

J,  Critical  J  toe  determiaiaf  tractate 

K/.Km  Mode-I  and  Mode-I  Street  iateaaity  (acton 

Ki ,  Critical  stress  iateaaity  factor  for  Mode-I  fractals 

A  Ratio  at  remote  lateral  load  to  remote  vertical  load 

r  Radial  diataace  from  crack-tip 

r,  Iacremeatal  crack  growth 

rc  Critical  crack  growth  iacremeat 

r,  Radios  deiaiag  core  regina 

5  Strata  eaergy  deaaity  factor 

S0,  S|,  S>  First  three  coefldeats  of  the  aeriea  expaaaioa  of  the  itraia  eaergy 

deaaity  factor 

S,  Strata  eaergy  deaaity  factor  at  load  iacremeat  j 

S,  Critical  itraia  eaergy  deaaity  facto*  for  determiaiag  fractare 

a  Iategratioa  path  for  the  J-integral 

T  Tractkw  vector  oa  a 

a  Displacement  vector  oa  a 

*i  Displacemeat  ka  the  y-directio#  at  tke  list  corner  node  behiad  the  crack-tip 

W  Strain  energy  density  faactioa 

W,  Critical  strain  eaergy  density  faactioa 

*,  |r  Cartesian  coordinate  axis  system 

s  Complex  vector  *  +  i  y 

a  Length  parameter  y/trfa 

9  Inclination  of  tke  crack 

r.r  Complex  faactioa  of  load,  biaxfahty  aad  angle  of  tar  bastion 

di  Diataace  of  the  list  comer  node  behind  the  crack-tip  to  the  crack-tip 

Normal  aad  abear  strains 
9  Angola*  coord iaate  at  crack-tip 

A  Faactioa  of  Poisson's  ratio 

¥  Pots  son's  ratio 

*  Ratio  of  the  drcamferesce  to  tke  d  is  meter  of  a  circle 

Noratal  aad  shea*  stresses 
Halomorphic  hartiaos 
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I.  INTRODUCTION 

The  aaatysis  of  the  itrm  held  aroaad  crack*  is  aa  elastic  solid  has  beea  wefl  doeaaseated  aad 
is  perhaps  nth  older  thaa  the  leid  sow  boss  as  Uaear  Elastic  Fraetare  Meehaaies.  Oae 
of  the  early  works  is  dae  to  Westerpaard  (1),  who  treated  the  sharp  crack  problem  aaiar  aa 
eifeaftuctioa  expaasioa  approach.  la  a  later  paper,  Weeterpaard  (3|  expressed  the  soiatioa  to 
the  saase  problem  la  terms  of  complex  aaatytk  fsaetkas.  This  later  work  has  siaee  fsiaed 
mack  recopiittoa,  aad  is  fteqaeatly  refsrtacid  by  the  aaalpsts  of Ractare  sssch safes  It  caa  be 
showa  that  the  Westerpaard  sslatioa  leads  to  the  followisf  classical  malt  lor  the  stress  held 
asar  the  tips  of  a  c  eat  re  crack  eoataiaed  ia  aa  iadaite  plate  loaded  by  a  remote  a  Biaxial  tea  dir 
strem  e  (see  Pip.  1): 


Ki  0 
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where  Ki  is  kaowa  as  the  strem  ist easily  (actor. 

However,  a  small  error  made  ia  the  Westerpaard  soiatioa  rraiaiaed  aadetected  lor  shanet 
thirty  years.  By  asiap  the  more  geaeral  complex  potest  iais  approach  of  MaskheiishviU  |S|, 
Sih  1 4)  showed  that  the  arbitrary  settiap  of  a  panic  alar  coastaat  ia  the  Westerpaard  soiatioa 
to  ssro  was  pearrafly  iavahd.  Discassioa  aa  the  effeeta  of  the  error  was  sabre qacatt>  takes 
ap  by  Eft  is  aad  Liebowiti  |S|,  aad  later,  EfUs  et  al  |«|  showed  thrn  the  mror  was  stahmlsal 
to  the  nmlmina  of  a  ana  slapitm  term  horn  the  streu  expremtnas.  By  eoasideilaf  aa  iadaite 
ceatre-c racked  plate  aader  biaxial  remote  loads  •  (perpeadiealar  to  the  crack)  sad  is  jpaaahrl 
to  the  crack),  it  was  showa  that  tke  soiatioa,  correct  to  the  moth  order  teem  ia  r.  Is 


(I*“) 
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Coaspariaf  Eqs  1.3  with  Eqs  1.1.  aad  settiap  i  -  0  lor  aa iaTial  loadiap.  H  is  aeea  that 
the  oaly  diBereace  is  the  a  aa -eta  palm  term  e  la  the  expeeaaioa  lor  »,  Of  coarse,  aa  r  —  0. 
this  teem  is  expected  to  become  asphfdde.  However,  it  don  sot  do  so  aa  rapidly  aa  oae  woaid 
like,  auksly  beesaar  of  the  aqaare  root  efect  of  the  tiapalar  term.  For  the  aaiaxtal  eaae  aa  aa 
example,  *  «  0  aad  Kt  ■  wjw i,  so  that  lor  •  ■  0*.  Eq.  1.3a  bee  oases 

11,1 

Em  the  error  eaaaed  by  asphetlaf  the  anaaiapsiaf  term  to  be  less  thaa  1%  say.  it  is  req aired 
that 


(Ml 
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-  SIxlO-4. 


(»•*» 


eHsascrilstiaa,lndstott*aEeaa Idhcrtpaa  das  with  ax*  sitrantst  data.  baaakesqantpaper, 
Mh «< alp|  <— oaatsatsd  tkat  the  local elastic  strab  imp  (hwtp  —4  at wh  wag  rate 
also  depead  ti*aldeiatly  oa  the  MaMpdlW  «M  load,  «4  tka  achhtary  oashaiaa  of 

Ik* - iafalac  um  Ml  tka  NkUM  to  tka  bMM  afleeta  nd  therefoae  leads  to  beasreet 

paodktMs  of  tka  ao«la  a f  growth  whea  asb*  Sik's  strata  aaartjr  dusky  eritaiioa  (•). 

Tfcts  anil  i  aattaaas  Iks  besbbbfaa  oa  tin  ifirl  irf  tfci  anlatiea  at  tka  — t  -Jir~*~r  * — 
ia  tka  fiimslatlna  of  tka  5  tkaary,  aad  tka  eaarr  genres  aad  dHnllks  la  Its  sppkeatioa  as  a 
steak.  It  vH  be  tkowa  tkat  area  tka  tkapia  aakadal  rasalts  rasp  ka  sicailcaatly  aftseted,  aad 
tka  aaatratt  psaetkkaar  of  tka  5  tkeory  eaa  ka  raialad. 


S.  •OUfTKNf  FOE  THE  INCUNED  CEACK  PEOBLEM 

Ik  Wastms  tka  egret  of  tka  aadasba  of  tka  aaa-ab*alsr  tern,  we  aril  Ian  to  a  paoklan 
far  wfcfek  a  wad  kaaara  ■ nielli a  exists.  Csastlsr  tka  Malta  plate  eoataiaiaf  a  eaatral  crack 
af  Isagth  la  sad  saki  sated  to  biaxial  leads  •  aad  h*  as  shows  la  FI*.  S.  Mowb*  Madtbel- 
IshrU  m.  kt  fan  nap  Me  tear  litas  wfclck  sapsaant  tka  saMjaa  ara 


(*.i) 

in*)  ■  </(#>  ■  j(ep  ♦  +  if ,  («) 

where 

if. -f<i  (!•>) 

j(»  +  r')-5(U-a**)  +  k<I  +  aM)}.  (td) 

aad  the  stresses  *..*»■*.»  art  pws  bp 

(l»l 


After  seats  aiaaipafatfaa.  sad  mil  rift*  Ike  safaiioa  sa*faa  to  0  <  r/a  <  I,  k  asap  he 
skawa  tkat  assstss  sapaastsa  af  tka  sskaids  ■  *t»aa 


I 


*i  17,  Ka.9U^  8„88\ 

*  1/ " ?&  *(  *  *) 

+  #(l  -6)eos2d, 


whet*  A"/  aad  Ka  ut  the 


iMmitjr  CKton  gives  by 


(J.7«) 

(276) 

(*.7c) 


if/  =  ^((i  +  *)-(t-*)««2/»|, 

jr,-^(l-*)riai^. 

Felmki  Hooke’s  taw,  the  ta-ptaae  elastic  itrmiat  are 

*»  ■  ]?(*»  “  ^fi)' 

«»  -  j(*»  -  A**), 

**»  JT 

f  «J  A  « «/  for  ptaae  stress, 

f  *  j-^j,  A  ■  j-^p  for  ptaae  strata, 

aad  E  aad  if  are  the  elastic  modstas  aad  Poisaoa's  ratio  respectively. 

Heaee,  the  strata  eaergy  deashy  faactioa  IP  for  a  ptaae  problem  is  gives  by 

*P  -  *(»*«.  +  #,«,  +  2r„e„) 

*  2?  I#*  +  #*  +  +  ~  '*'»)] 

Sabsthatiag  Eqs  2.7  iato  Eq.  2.11,  we  obtaia 


(2.8) 

(2-9) 


(2.10) 


(2.11) 


IP  -  -  +  4=  +  C,  (2.12) 

r  y/f 

where 

4  “  wk<W?,1  +  eos#)()S -cos#)  -  All  +  cos 8)1 

+  X*((9  -  4 cos#  +  Scos28)  +  A(1  +  4eos#  4-  2 cos 2# )| 

+  «IT/X»staH(eos8-l)  +  A(eos*  +  l)|}, 

8  *  — ~  {AT(  cos  J((2  -  cos#  +  coslf )  -  A(2  +  cos#  -  eos2f)| 

If  v*t  • 

-  Ka  sta  J|(4  +  eo»8  +  cot  28)  +  A(coe8  +  cos 28))}, 

C  -  ^^(l-fcj’eos^i. 
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Not*  that  tht  iacbaio*  at  the  »oa  *pla  tana  la  th*  Mm  ray  arnica  bads  to  two 
additional  Una  hi  IT,  at*  at  whkh  ha  a  r~lf*  kaqabri*  whfct  th*  othat  It  a»M(ikr. 
TO  mo  how  thk  dheti  th*  appMcatba  of  th*  f  theory,  It  i«  batractiv*  to  aanniriap  th*  bosk 
of  th*  thacay. 


2.1  Th*  Strata  foirp  Polity  Crftorlow 

Th*  strata  energy  danaity  criterion  wm  tnt  proposed  by  Sih  (•),  aad  ita  appfcatio*  ha*  beta 
w*S  docamawted  (see  to*  eaafb  |9j  -  |11|).  It*  fseantaHs*  b  baaed  oa  th*  aoaaa^tlo*  that 
tike  tttahi  «a**gy  dtaaity  faaetioa  aear  the  crack-tip  possesses  a  1/r  atagabrity.  Thk  aeay  be 
•eea  hoot  Eq.  2.12  ia  whkh  for  saHeieatly  total  r,  W  stay  b*  written  a* 

*  -  7.  (*•»*) 

■•ace,  aa  (-independent  parameter,  kaowa  at  th*  tttaia  eaetgy  deaaity  (actor  5,  May  be 
detaed  each  that 


S*rff  *  A.  (2.14) 

The  ciiterioa  concerning  crack  iaitiatkoa  aad  the  aagb  at  growth  k  thea  baaed  oa  the 
iolowiag  three  hypotheses: 

Hypotheak  X:  Crack  extension  begins  aloag  the  direction  where  the  strata  energy  density  (actor 
k  a  local  tntahwam  (deaoted  by 

Hypotheak  2:  Fractare  it  homiaeat  whea  the  local  atiahaan  Smtm  (cache*  a  critical  valae  S., 
aad  that  5,  it  a  material  parameter  which  characterises  the  (ractnre  atreagth 
of  the  material. 

Hypotheak  S:  The  aamaat  at  stabh  iacremeatal  growth  ri.r*,.. ,r,  la  governed  by 

contl.  (2.15) 

*i  *a  *j  r,  '  ’ 

where  a  poo  reach  lag  r„  aaatabie  fractare  ocean. 

For  a  ceot re-cracked  apeebaea  oader  mode  I  eoaditioaa,  both  the  directioo  of  crack  prop- 
agatioo  aad  fractare  tooghaen  are  well  kaowa.  Heace,  the  aabatitatioB  of  t  »  0*,ff/  = 
Ki.,  Kg  mO  iato  Eq.  2.14  give*,  tor  plaae  strata  eoaditioaa,  a  coaveaieat  expression  for  relat- 
iag  S.  to  existing  fractare  toaghaeaa  data,  via., 


,  (l  +  i/)(l-2k)  , 
5*  - vTr - Ki. 


(*•!•) 


Figare  2  shows  a  graphical  repreaeatatioa  of  Hypotheak  2.  It  show*  that  the  amooat  at 
lacreaaeatal  growth  eaa  be  predicted  by  the  iatenectioa  of  the  plots  of  IF  aad  the  material 
coastawt  IF..  Ia  aa  ehatfe  aaalysk,  W  k  aaboaaded  a*  r  -»  0  for  aay  ghrea  applied  atreas  #  so 
that  IF  wH  ahrays  tatervect  IF..  la  practice  however,  a  crack  wlB  aot  grow  oaten  a  saHcieat 
load  k  applied.  1b  aecoaat  for  thk  phenomena*,  a  threshold  distaaee  r,  k  deiaed  sack  that 
fce  r  <  e.  (referred  to  as  the  core  regioa),  coatlaaam  meckaaics  k  assnmed  to  (ail  short  of 
reality  so  that  aay  aaalytk,  aad  heace  predict ioa  of  growth,  aaest  be  hep*  oataide  the  cose 
regioa.  That  k,  the  coadkioa  rt  >  r,  atast  be  satkSed.  It  k  abo  dear  that  tor  a  material 
whkh  exhibits  aegHgibb  atabb  growth  (Le.,  a  brittle  ataterial),  r.  «  re.  For  a  4140  steel  with 
variooa  yield  atraagtha  aad  fractare  toaghaesaea,  Sik  |1H|  showed  that  r«  raaged  from  0.00005 
ia.  (0.0146  mm)  to  0.01246  ia.  (0.2416  mm). 
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3.  THE  APPLICATION  OF  THE  S-THEORY  IN  PREDICTING  BRITTLE  FAIL¬ 
URE 

One  advantage  of  the  strain  energy  density  criterion  over  most  other  fracture  criteria  is  that 
both  the  angle  of  growth  and  critical  load  may  be  predicted  by  a  single  parameter.  Furthermore, 
without  the  restriction  of  seif-similar  growth  as  required  by  other  criteria  such  as  the  critical 
energy  release  rate  approach,  the  5  criterion  may  be  applied  in  complex  load  systems.  For 
example,  by  means  of  finite  element  modelling,  the  5  distribution  around  a  crack-tip  contained 
in  a  complex  structure  under  an  arbitrary  load  system  may,  in  theory,  be  used  to  predict  the 
failure  load  as  well  as  the  direction  of  propagation.  However,  it  will  be  shown  in  this  section 
that  without  sufficient  care,  this  practice  can  lead  to  erroneous  results. 


3.1  Predicting  the  Onset  of  Failure 

Consider  the  uniaxial  problem  described  previously.  Since  the  direction  of  growth  for  this  case 
is  known  a  priori  (namely  along  the  x-axis),  substituting  9  =  0“,  ff  =  90*,  Ku  =  0  into  Eq.  2.12 
and  multiplying  through  by  r  yields  an  expression  for  5  along  the  critical  direction,  viz., 


^(1-  A)  -X)J 

2  r  2 

—  +<r2r 

x  V 

X 

(3.1) 


And  since  K/  —  o^/ia. 


S  =  ^Ml-A)-(l-A)v/2^  +  r] 


(3.2) 


The  one-term  representation,  Eq.  2.14,  may  be  obtained  by  taking  the  limit  of  Eq.  3.2  as 
r/a  —  0,  giving 


S  = 


»“-*>■ 


(S-3) 


Taking  a  typical  value  of  v  =  0.3,  Fig.  4  shows  a  comparison  between  the  one-term  rep¬ 
resentation  and  the  higher  order  expression  for  5  .  It  may  be  seen  that  there  is  significant 
difference  between  Eq.  3.2  and  Eq.  3.3.  At  the  relatively  small  value  of  r/a  =  0.02,  the  neglect 
of  the  higher  order  terms  gives  rise  to  errors  of  approximately  20%  in  both  the  plane  stress 
and  plane  strain  results.  Indeed,  values  of  r/a  >  0.02  have  not  been  uncommon  in  many  anal¬ 
yses  using  the  strain  energy  density  method  (e.g.,  |10),  (13)).  What  is  more  important  in  this 
example  is  that  along  the  critical  direction,  the  one-term  expression  always  over-estimates  the 
correct  valne  of  5  for  r/a  >  0.  Since  5,  is  derived  from  the  one-term  expression  (Eq.  2.16), 
the  prediction  of  fracture  load  using  the  5-theory  by  means  of  finite  element  modelling,  or  by 
physically  monitoring  the  elastic  strains  at  some  small  but  flnite  distance  from  the  crack  tip, 
will  invariably  produce  an  over-estimate.  In  the  above  example,  if  the  strain-energy  density  5 
at  r/a  =  0.02  is  monitored,  and  Hypothesis  2  of  the  criterion  is  applied  so  that  fracture  load 
is  assumed  to  be  reached  when  S  =  5e,  then  S(r/o  -»  0)  would  be  greater  than  5C  by  about 
20%  and  hence  the  fracture  strength  would  be  over-estimated  by  9.5%.  The  non-conservative 
nature  and  the  relatively  large  magnitude  of  the  error  involved  is  certainly  undesirable,  and  it 
shows  how  the  blind  application  of  the  5-theory  in  design  work  can  be  extremely  dangerous. 

For  the  5-theory  to  be  applicable,  the  analysis  must  be  confined  to  a  region  for  which  the 
one-term  expression,  Eq.  3.3,  is  valid.  Allowing  an  error  of  say  5  %,  we  get  from  Eq.  3.2 


<  0.05. 


(J.4) 


(1  -  X)<J2r/a  -  r/a 
(l-X)-(i-X)^/2rJa  +  r/a 

The  solution  which  satisfies  the  coeditiom  0  <  r/a  <  1  it  ghreu  by 

T-  <  0.00124.  (3.5) 

Equation  3.5  represents  am  upper  limit  to  the  region  where  the  strain  energy  density  analysis 
may  be  applied  with  acceptable  accuracy.  However,  we  recall  that  there  exists  also  a  lower  limit 
r„  which  defines  the  core  region  (see  $2.1).  Consequently,  the  region  where  the  5-theory  is  valid 
may  be  given  by 


r„  <  r  <  0.00124  a. 


(3.6) 


Equation  3.6  implies 


“  “  0.00124' 


(3.7) 


The  restrictions  inferred  by  Eqs  3.5  and  3.7  are  severe,  and  in  some  cases,  may  be  considered 
as  impractical.  For  the  high  fracture  strength  4140  steel  ( r0  =  0.01345  in.)  quoted  in  |12)  as  an 
example,  the  minimum  crack  length  required  for  a  valid  LEFM  analysis  would  be  10.85  in.  (275.6 
mm).  Such  a  condition  would  be  difficult  to  satisfy,  particularly  for  real  life  situations  where 
cracks  of  much  shorter  lengths  are  detected  and  require  analysis.  Equation  3.5  on  the  other 
hand  can  almost  always  be  satisfied  in  a  finite  element  model  in  which  mesh  sizes  of  the  order 
of  0.001  a  may  be  readily  handled  on  a  modern  computer.  However,  using  such  a  fine  mesh 
for  a  linear  analysis  may  be  difficult  to  justify  as  stress  intensity  factors  have  been  successfully 
computed  by  other  methods  with  much  coarser  grids.  To  highlight  this  point,  a  finite  element 
model  of  a  centre-cracked  plate  (i/  —  0.3)  subjected  to  uniaxial  tensile  loading  and  plane  strain 
conditions  was  established.  Because  of  symmetry,  only  one  quarter  of  the  plate  which  has  a 
width  to  height  ratio  6/A  =  1,  and  a  crack-length  to  width  ratio  a/6  =  0.25  was  modelled. 
Figure  5  shows  the  two  mesh  schemes  adopted  for  the  analysis.  The  parameters  calculated 
include 

1)  The  stress  intensity  factor  given  by  the  expression  (14] 


where  vt  and  Ai  are  respectively  the  y-displacement  and  distance  from  the  crack  tip  at 
the  first  corner  node  behind  the  crack  tip,  and  F  is  as  defined  in  Eq.  2.10. 

2)  The  7-integral  given  by  the  numerical  integration  of  the  expression  (15) 

J  =  Jwdy-T-^d',  (3.9) 

where  T  and  u  are  the  traction  force  and  displacement  vectors  along  the  path  t  respectively. 
The  integration  path  chosen  in  this  case  is  shown  in  Fig.  5.  It  may  be  shown  that  7  is 
path-independent  for  an  elastic  material  and  b  related  to  the  stress  intensity  factor  by  the 
following  expression 


3)  The  strain  energy  density  factor  S(=  rW)  at  nodal  points  along  the  critical  direction 
(3  =  OP),  where  W  is  calc  slated  from  the  computed  stresses  (Eq.  3.11).  To  estimate  the 
stress  intensity  factor  from  the  strain  energy  density  factor,  the  conventional  expression 
which  neglects  the  higher  order  terms  is  used  (12) 


|,u| 

In  the  following  presentation,  all  values  of  K/  and  5  have  been  normalised  by  Oy/a  and 
a* aft  respectively.  The  computed  strain  energy  density  along  the  critical  direction  (•  =  0*) 
for  Mesh-1  is  shown  in  Fig.  6.  An  immediate  observation  is  the  systematic  scatter  of  the  results 
about  some  mean  valae.  The  strain  energy  density  factor  appears  to  be  overestimated  at  the 
corner  nodes,  and  under-estimated  at  mid-side  nodes.  Whilst  the  average  deviation  of  5  from 
the  mean  amounts  to  appradmately  S.7  %,  the  scatter  of  the  computed  stresses  would  be  only 
approximately  2.8  %  and  is  therefore  considered  acceptable.  If  the  assumption  of  a  constant  S  is 
made  over  the  region  say  0  <  r/a  <  0.04,  then  a  best  fit  to  the  data  gives  the  result  S  -  0.2597, 
which  when  used  in  Eq.  3.11,  gives  K /  =  1.771.  This  compares  with  Ki  =  1.950  obtained  using 
the  analytical  solution  given  in  Rooke  and  Cartwright  [16],  and  therefore  represents  a  9.2  % 
underestimation. 

The  computational  results  for  the  reined  mesh  are  shown  in  Fig.  7.  Again,  a  similar 
systematic  scatter  is  apparent,  and  the  assumption  of  a  constant  5  in  the  range  0  <  r/a  <  0.04 
gives  5  =  0.2837,  or  from  Eq.  3.11,  Ki  =  1.848  which  is  in  error  by  5.2  %.  However,  the 
downward  trend  at  the  5-distribution  as  predicted  by  the  analytical  solution  is  clearly  evident 
in  this  case,  and  compares  well  with  the  plot  of  Eq.  3.1  (where  If  j  is  taken  to  be  1.950).  This 
suggests  that  some  sort  of  curve  Itting  to  the  correct  form 

S  =  S.  +  S,^  +  SiQ,  (312) 

would  provide  a  more  accurate  prediction  of  strain  energy  singularity.  This  was  indeed  found 
to  be  the  ease.  Applying  a  least  squares  It  for  the  5  data  of  Mesh-2  to  Eq.  3.12,  it  was  found 
that  So  =  0.3082.  Substituting  Sa  for  S  in  Eq.  3.11  yields  Ki  =  1.930,  giving  an  error  of  only 
1%. 

Table  1  summarises  the  various  methods  of  computing  the  stress  intensity  factor  and  their 
errors  with  respect  to  the  analytical  solution  obtained  in  (16). 


K, 

Ret  |l«j 

Ki(v„A,) 

K,(J) 

Jf;(5  =  court.) 

Mesh-2 

ia 

1.892 

_ 

u 

l 

1.894 

i 

1.896 

1.771 

i 

1.848 

%  Error 

2.5 

ma 

2.8 

5.2 

Table  l.  Stress  Intensity  Factors  Calculated  by  Various  Methods 


7 


It  may  be  seen  from  table  I  that  the  stress  intensity  {actor  is  predicted  accurately  by 
both  the  crack  opening  and  the  /-integral  even  for  Mesh-l.  The  evaluation  of  A'/  from  5 
however  required  the  much  finer  mesh  to  achieve  an  accuracy  to  within  5  %.  For  a  more 
accurate  prediction,  a  least  squares  fit  of  the  correct  form  was  required.  Whilst  this  curve 
fitting  technique  proved  to  be  useful  in  the  case  of  Mesh-2,  it  would  be  unreasonable  to  apply 
this  to  the  data  of  Mesh-1  where  a  poor  correlation  to  the  true  solution  is  expected. 

It  is  also  interesting  to  point  out  that  whilst  the  /-integral  involves  the  integration  of  a 
strain  energy  term,  the  problems  associated  with  the  5-theory  do  not  appear,  as  it  has  been 
shown  by  Eftis  et  aJ  [7|  that  the  truncation  of  the  higher  order  terms  has  no  effect  on  /  and 
Eq.  3.10.  Another  feature  of  /  is  that  the  integration  process  tends  to  nullify  the  systematic 
scatter  of  the  numerical  data  and  was  therefore  able  to  predict  Ki  accurately  without  the 
application  of  any  best  fit  technique.  In  practice,  it  is  most  likely  that  /  is  used  in  a  non¬ 
linear  analysis  where  stress  intensity  factors  are  not  defined.  Hence,  /  is  not  generally  used  for 
computing  Ki,  but  instead,  it  is  calculated  and  compared  directly  to  a  critical  value  Jc  (which, 
like  Kit  and  Sc,  is  a  material  parameter)  to  determine  whether  failure  is  to  occur.  The  above 
results  merely  show  that  a  finite  elements  approach  to  the  /-integral  can  be  applied  with  some 
confidence.  However,  it  must  be  remembered  that  the  restriction  of  self-simitar  growth  must  be 
observed,  and  that  its  usefulness  when  extended  to  true  elasto-plastic  materials  is  yet  uncertain 
|17|. 


S.I  Predicting  the  Direction  of  Growth 

Having  seen  how  the  neglected  higher  order  term  can  affect  the  prediction  of  crack  initiation  by 
the  5-theory,  its  effect  on  the  prediction  of  the  direction  of  crack  extension  is  now  examined.  In 
testing  the  maximum  normal  tensile  stress  theory  for  predicting  the  direction  of  crack  extension, 
Williams  and  Ewing  (18]  conducted  experiments  on  PMMA  (polymethylmethacrylate)  sheet 
specimens  which  contain  a  central  inclined  slit  crack  and  subjected  to  uniaxial  tensile  loads. 
It  was  found,  particularly  for  steep  crack  angles  (/)  — >  0*),  that  the  angle  of  growth  deviated 
from  the  predicted  results  presented  by  Erdogan  and  Sih  [19|.  The  authors  attributed  the 
discrepancies  to  the  fact  that  only  a  one-term  expansion  was  used  in  the  expression  for  the 
stresses,  and  showed  that  by  selecting  a  critical  parameter  a  =  \j2r/<i  =  0.1(r/a  =  0.005),  and 
including  higher  order  terms  in  the  expansion,  a  better  fit  of  experimental  data  with  theory  is 
achieved.  In  a  subsequent  discussion,  Erdogan  and  Sih  |20]  showed  that,  unlike  the  maximum 
normal  tensile  stress  criterion,  the  strain  energy  density  criterion  provides  a  better  fit  to  the 
experimental  data  and  is  relatively  insensitive  to  the  parameter  or.  To  investigate  this  point 
further,  Eq.  2.12  is  differentiated  with  respect  to  9  to  form 


as  _  dA  aB  r 

a^s  =  aM  a3a  r 
ae3  ~  a«s  +  a«s  v/r' 


(3.13) 


and  from  Hypothesis  1  of  the  criterion,  the  direction  of  crack  propagation  is  determined  when 


and 


(3.14) 


Assuming  u  =  0.33  and  for  plane  strain  conditions,  the  solutions  to  Eq.  3.14  over  the  range 
0*  <  9  <  90*  for  a  =  0,0.1  and  0.2  are  presented  in  Fig.  8  together  with  the  scatter  band  of  the 
results  in  [IS].  It  may  be  seen  that  the  one-terra  representation  of  the  S-theory  (equivalent  to 
the  case  or  =  0)  fitted  the  experimental  data  better  than  the  one-term  approach  to  the  maximum 
stress  criterion.  It  is  also  noted  that  the  solutions  for  both  a  =  0.1  and  0.2  fall  within  or  near  the 
experimental  scatter  band.  However,  there  is  some  evidence  that  the  use  of  a  =  0  for  0  <  45* 
and  a  —  0.2  for  0  >  45*  can  result  in  a  more  accurate  prediction.  This  has  the  implication  that 
the  core  radius  r„  may  in  fact  not  be  a  material  property.  The  question  of  whether  r„  may  be 
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considered  as  a  material  constant  has  previously  been  addressed  by  Chang  {31  j.  Chang  found 
that  in  order  to  obtain  reasonable  agreement  between  the  experimental  data  o it  |18),  |22|,  |23| 
and  the  strain  energy  density  theory,  various  values  of  rc/a  (ranging  from  0.005  to  0.15)  had  to 
be  used  despite  the  fact  that  all  data  were  presented  for  the  same  material  (PMMA)  and  for 
similarly  sized  cracks.  He  therefore  concluded  that  r0/a  (and  hence  r0)  ‘can  hardly  be  justified 
as  a  material  parameter  in  the  5-theory'. 

Chang  [21]  also  discussed,  at  length,  the  dilemmas  which  may  arise  in  applying  the  5- 
theory  when  no  relative  minimum,  or  alternately,  when  more  than  one  relative  minimum  in  5 
exist  in  the  solution.  Swedlow  [24]  showed  that  tor  uniaxial  loading  configurations,  the  choice 
of  the  global  minimum  leads  to  incorrect  predictions.  Swedlow  then  proposed  the  additional 
requirement  that  the  5m,„  which  governs  fracture  must  be  associated  with  a  tensile  hoop  stress. 
On  the  other  hand,  Sih  and  Madenei  [10]  assert  that  it  is  the  maximum  of  all  5m,,’r  which 
first  reaches  the  material  threshold  and  is  therefore  the  critical  factor.  For  the  inclined  crack 
problem  considered,  two  local  minima  are  found  for  all  0  except  for  0  —  90*,  and  indeed,  it  is 
the  maximum  of  the  two  at  any  given  r  which  corresponds  to  the  results  presented  in  Fig.  8. 
However,  an  interesting  situation  can  arise  when  other  loading  conditions  are  considered.  As  an 
example,  under  a  biaxial  tension-compression  loading  system  with  k  =  -1,  and  for  0  =  60*,  the 
paths  of  the  two  minima  are  as  shown  in  Fig.  9.  The  corresponding  plots  of  S  along  the  paths 
denoted  by  i  and  ii  are  presented  in  Fig.  10.  It  is  clear  from  the  plots  that  the  determination 
of  the  maximum  of  the  minima  would  depend  on  the  selection  of  r/a. 

It  has  been  seen  in  the  uniaxial  case  that  although  the  use  of  the  5-theory  to  predict 
crack  initiation  requires  an  extremely  small  value  of  r/a,  the  restriction  was  much  less  severe 
for  predicting  the  propagation  angle  as  reasonable  predictions  may  be  achieved  for  the  range 
0  <  r/a  <  0.02.  However,  the  above  example  shows  that  this  may  not  be  valid  in  general  as  the 
choice  of  r/a  appears  to  be  crucial  in  determining  whether  the  crack  is  to  grow  in  direction  i  or 
ii.  Perhaps  this  is  also  evidence  for  the  possibility  that  the  5-theory  alone  may  be  insufficient 
in  determining  the  direction  of  crack  propagation  in  general,  and  a  modification  such  as  that 
proposed  by  Swedlow  [24]  may  be  in  order. 


4.  CONCLUSION 

It  has  been  shown  that  the  usual  assumption  of  a  1/r  energy  singularity  is  valid  only  within 
an  extremely  small  regime  around  the  crack  tip.  For  the  centre-cracked  plate  considered,  this 
region  is  typically  of  the  order  of  r/a  <  10~s.  As  a  consequence,  the  application  of  the 
5-theory  at  some  distance  outside  this  region  may  result  in  substantial  errors  in  the  prediction 
of  crack  initiation.  For  a  finite  element  analysis,  an  extrapolation  technique  using  the  more 
accurate  form  is  proposed  and  found  to  be  useful.  However,  this  technique  still  requires  a 
relatively  fine  computational  grid,  and  a  guideline  for  maximum  mesh  size,  optimum  number  of 
data  points  and  the  valid  extrapolation  domain  for  arbitrary  crack  and  loading  configurations 
has  yet  to  be  established. 

It  has  also  been  shown  that  the  restriction  on  r/a  is  somewhat  less  severe  for  predicting 
the  crack  growth  direction,  and  that  for  the  uniaxial  load  case,  reasonable  agreement  between 
experimental  data  and  predicted  results  may  be  achieved  over  a  relatively  large  range  of  r/a. 
Unfortunately,  this  may  not  be  taken  as  a  general  rule  as  illustrated  by  the  biaxial  load  example 
where  a  dilemma  iu  choosing  the  correct  5m(,  may  arise  when  r/a  is  arbitrarily  selected. 

In  closure,  it  should  be  emphasised  that,  despite  the  problems  revealed  by  the  current 
work,  the  strain  energy  density  factor  should  not  be  disregarded  as  a  useful  parameter.  There 
is  no  doubt  that  the  5-theory  works  well  under  certain  conditions,  and  its  potential  in  handling 
mixed  mode  fractures  is  particularly  valuable.  The  close  relationship  between  S  and  the  stress 
intensity  factors  (Eq.  2.12),  and  the  reasonable  agreement  between  predicted  and  existing 
experimental  data  on  propagation  angles,  teud  to  support  this.  On  the  other  hand,  limitations 
to  the  theory  must  be  identified  and  realised.  Questions  such  as  whether  or  not  r0  is  a  valid 
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material  parameter,  or  what  shoaid  be  doae  whew  mahipie  miaima  ia  5  exist,  have,  ia  tbe 
preteat  aatbor'i  optaioa,  yet  to  be  positively  resolved. 
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Until  recently,  the  one-parameter  singular  expression  for  stresses  near  a 
crack-tip  was  widely  thought  to  be  sufficiently  accurate  over  a  reasonable  region  for 
any  geometry  and  loading  conditions.  This  view  has  been  fast  changing  due  to  the 
growing  evidence  that  the  inclusion  of  higher  order  terms  can  significantly  affect  the 
solution,  particularly  under  certain  biaxial  loading  conditions.  In  this  context,  the 
present  paper  examines  the  strain  energy  density  criterion  for  fracture,  and  the 
consequences  of  the  assumption  of  a  1/r  energy  singularity  in  the  formulation  on  its 
application.  It  is  found  that  this  assumption  imposes  a  rather  severe  restriction  on  the 
region  for  which  the  criterion  is  applicable,  and  that  its  application  on  an  arbitrarily 
selected  ’small’  distance  from  the  crack-tip  (a  procedure  which  has  been  adopted  by 
many  experimentalists  and  finite  element  analysts),  can  lead  to  erroneous  results. 
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